Film-thickness dependence of 10 GHz Nb coplanar-waveguide resonators 
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We have studied Nb A/2 coplanar-waveguide (CPW) resonators whose resonant frequencies are 
10 — 11 GHz. The resonators have different film thicknesses, t = 0.05, 0.1, 0.2, and 0.3 /xm. We 
measured at low temperatures, T = 0.02 — 5 K, one of the scattering-matrix element, S21, which is 
the transmission coefficient from one port to the other. At the base temperatures, T — 0.02 — 0.03 K, 
the resonators are overcoupled to the input/output microwave lines, and the loaded quality factors 
are on the order of 10 3 . The resonant frequency has a considerably larger film-thickness dependence 
compared to the predictions by circuit simulators which calculate the inductance of CPW taking 
into account L g only, where L g is the usual magnetic inductance determined by the CPW geometry. 
By fitting a theoretical S21 vs. frequency curve to the experimental data, we determined for each 
film thickness, the phase velocity of the CPW with an accuracy better than 0.1%. The large 
film-thickness dependence must be due to the kinetic inductance Lk of the CPW center conductor. 
We also measured S21 as a function of temperature up to T — 4 — 5 K, and confirmed that both 
thickness and temperature dependence are consistent with the theoretical prediction for Lk- 
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I. INTRODUCTION 

Microwave resonators (for example, Chap. 7 of Ref. [l|) 
are one of the key components in a variety of circuits 
operated at GHz frequencies, and their new applications 
continue to emerge. A simple example is band-pass fil- 
ters, which are based on the fact that the microwave 
transmission through resonators is frequency sensitive. 
The same idea is also used for more complex devices, 
such as oscillators, tuned amplifiers, and frequency me- 
ters. Actually, having high-quality filters and oscilla- 
tors is critical in mobile communications, where avail- 
able bands keep getting overcrowded as demand grows 
rapidly. 

Another application of microwave resonators is radia- 
tion detectors, which often consist of sensor heads and 
read-out circuits, and resonators can be used in the 
readout circuit. When one would like to detect at the 
single-photon level, one needs to have detectors with 
high enough energy resolutions. In this respect, super- 
conducting sensor heads^i 3 -^ can be advantageous, and 
may be the only solution at present depending on the 
energy range of the object. Once one decides to use su- 
perconducting sensor heads, it makes sense to fabricate 
the read-out circuit with superconducting materials as 
well. Superconducting microwave resonators allow one 
to obtain higher quality factors, which are favorable for 
frequency multiplexing. In addition, one type of photon 
detector is designed to probe the change in the kinetic in- 
ductance of superconducting thin-film resonator due to 
the absorbed photons^ In this device concept, the res- 
onator works as a sensor head rather than a part of the 
readout circuit. 

Recently, superconducting resonators are used for 
the nondemolition readout of superconducting qubits as 



welli^ Since the demonstration by Wallraff et al.,® this 
type of readout scheme has been one of the main topics in 
the field of superconducting qubits, and we are also devel- 
oping a similar readout technique^ In superconducting 
resonators, kinetic inductance, which is essentially the 
internal mass of the current carriers, plays an important 
role especially when the superconducting film is thin. In 
our circuit j£ for example, a Nb A/4 coplanar-waveguide 
(CPW) resonator is terminated by an Al dc SQUID, and 
the total thickness of the Al layers is 0.04 fim. In or- 
der to avoid a discontinuity at the Al/Nb interface, we 
usually choose the Nb thickness to be 0.05 /zm, which is 
much thinner than a typical thickness of > 0.3 /jm for su- 
perconducting integrated circuits fabricated by the stan- 
dard photolithographic technology. Fabricating circuits 
with thinner films is actually important from the view- 
point of miniaturization as well. Therefore, for designing 
resonators, quantitative understanding of the kinetic in- 
ductance in the CPW is important. 

There have been a number of reports on kinetic in- 
ductance for a variety of materials ^' 9 ' 10 ! 11 ] 12 ? 13 In gen- 
eral, however, kinetic inductance is indirectly measured 
by assuming a theoretical model, and as a result, the 
uncertainties are relatively large. Thus, although ki- 
netic inductance is a well established notion and the phe- 
nomenon is qualitatively understood, the quantitative 
information is not necessarily sufficient from the point 
of view of applications, especially at high frequencies, 
>10 GHz. When we would like to precisely predict the 
resonant frequency, the best solution would be to char- 
acterize the actual CPW in a simple circuit. Such char- 
acterization should also improve the knowledge of super- 
conducting microwave circuits. Very recently, Goppl et 
aliU measured a series of Al CPW resonators with nomi- 
nally the same film thickness of 0.2 /xm, and investigated 
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TABLE I: List of resonators, t is the thickness of Nb film; f r 
is the resonant frequency; Ql is the unloaded quality factor; 
C c is the coupling capacitance; v p is the phase velocity, and 
its ratio to the speed of light c is listed in percent. The values 
for f r and Q l are obtained at the base temperatures. C c and 
v p are evaluated by least-squares fitting (see Fig. [4} with C = 
f.6 x 10~ 10 F/m, and their uncertainties are determined by 
changing the value of C by ±10%, where C is the capacitance 
per unit length. 
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FIG. 1: Schematic diagram of coplanar-waveguide (CPW) 
resonators. A CPW of length / is coupled to the microwave 
lines through capacitors C c . 



the relationship between the loaded quality factor at the 
base temperature of 0.02 K and the coupling capacitance. 
For this purpose, it is justified to neglect kinetic induc- 
tance because the kinetic inductance should be the same 
in their resonators and estimated^ to be about two or- 
ders of magnitude smaller than the usual magnetic induc- 
tance determined by the CPW geometry. In this work, 
on the other hand, we paid close attention to the resonant 
frequency as well, and characterized Nb CPW resonators 
as a function of film thickness rather than a function of 
coupling capacitance. We also looked at the tempera- 
ture dependence in order to discuss kinetic inductance in 
detail. 



II. EXPERIMENT 

We studied two series of Nb A/2 CPW resonators listed 
in Table HI Each resonator consists of a section of CPW 
and coupling capacitors, as shown schematically in Fig.[T] 
The resonators were fabricated on a nominally undoped 
Si wafer whose surface had been thermally oxidized. On 
the SiC>2/Si substrate, a Nb film was deposited by sput- 
tering and then patterned by photolithography and SFq 
reactive ion etching. Figure HJa) represents the cross 
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FIG. 2: (Color online) (a) Cross section of a coplanar waveg 
uide. (b) Top view of coupling capacitors. 
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FIG. 3: Schematic diagram of a typical measurement setup. 
Boxes represent attenuators. 



section of CPW. The center conductor has a width of 
w = 10 /im, and separated from the the ground planes 
by s = 5.8 /im, so that the characteristic impedance be- 
comes - 50 n. The thickness of Nb is t = 0.05, 0.1, 0.2, 
or 0.3 /im (see Table QJ, an d that of SiC^/Si substrate is 
h = 300 /im. The Si02 layer, whose thickness is 0.3 /im, 
is not drawn in Fig. [2ja). We employed interdigital cou- 
pling capacitors as shown in Fig.^b). The finger width 
is Wf = 9 /im, the space between the fingers is s/ = 2 /im, 
and the finger length is If = 78 /im for Resonators A1-A4 
and = 38 /im for Resonators B1-B4. Here, we quoted 
designed dimensions for the Nb structures. The actual 
dimensions differ by about 0.2 /im due to over-etching; 
for example, w and Wf are ^0.2 /im smaller, whereas s 
and Sf are ~ 0.2 /im larger. In this paper, we define the 
resonator length I as the distance between the center of 
the fingers on one side and that on the other side, and 
I = 5.8 mm for all resonators. Because our chip size is 
2.5 mm by 5.0 mm, our CPWs meander as in Fig. [3] 

The resonators were measured in a 3 He- He dilution re- 
frigerator at T = 0.02 — 5 K. A typical measurement setup 
is shown schematically in Fig. [3] The boxes in the figure 
represent attenuators. The amount of attenuation was 
not the same because the microwave lines in our refrig- 
erator had been designed for several different purposes. 
The attenuation was x = 10 dB for Resonators A2, B3, 
and B4, and x = 20 dB for the others; y — 10 dB for all 
resonators except Al and A3. For Resonators Al and A3, 
we used a line with no attenuators (y = dB) but with 
an isolator and a cryogenic amplifier at 4.2 K. The gain of 
the cryogenic amplifier was 40 dB for Resonator Al and 
34 dB for A3. We measured the transmission coefficient 
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FIG. 4: (Color online) Amplitude of the transmission coeffi- 
cient S21 as a function of frequency for (a) Resonators A1-A4, 
and (b) Resonators B1-B4. 



£21 by connecting a vector network analyzer to the "IN" 
and "OUT" ports in Fig. [3] A typical incident power 
to the resonator was —40 dBm. For each resonator, we 
confirmed that the measurements were done in an appro- 
priate power range in the sense that the results looked 
power independent. 



III. RESULTS 
A. S21 at the base temperatures 

Figure 0] shows the amplitude of 52 1 at the base tem- 
peratures, T = 0.02 — 0.03 K, as a function of frequency 
/ for all resonators. The resonant frequency f r has a 
rather large film-thickness dependence. Our interpreta- 
tion is that this is due to the kinetic inductance of the 
CPW center conductor. Before discussing the thickness 
dependence in detail, let us look at the quality factors. 

What we obtain by measuring 52 1 as a function of / 
is the loaded quality factor Ql, which is related to the 
external quality factor Q e and the unloaded quality factor 
Qby 



Qi 



Qe+Q- 



(1) 



In general, Q e is determined mainly by C c , whereas Q is 
a measure of the internal loss, which arises not only from 
the dielectric but also from the superconductor in the 
high-frequency regime. Our resonators should be highly 
overcoupled to the input/output lines at the base tem- 
peratures, that is, Q » Q e , and thus, Ql ~ Qe- As 
listed in Table HI Ql of our resonators is on the order of 
10 3 . These values are not only reasonable for the designs 
of our finger-shaped coupling capacitors but also much 
smaller than typical values of Q below 0.1 K for super- 
conducting microwave resonators ^ 13 ' 14 When Q > Q e , 
the maximum |52i| is expected to be dB. We have 
confirmed by taking into account attenuators, amplifiers, 
and cable losses, that our measurements are indeed con- 
sistent within the uncertainties of gain/loss calculations, 
1-2 dB. Based on this confirmation, the experimental 
data in Fig. 0] are normalized so that the peak heights 
equal dB. 



The solid curves in Fig. [4] are calculations based on the 
transmission (ABCD) matrix (for example, Sec. 5.5 of 
Ref. [l|), and they reproduce the experimental data well. 
The matrix for the resonators is given by 
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j is the imaginary unit, 



cos (31 



jZ cpw sin (31 
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(4) 



lcpw ~ V j{Z epw )- x ^D.pl cos (31 
for lossless CPWs, to = 2nf, (3 — ui/v p , 

v p = 1/VZC (5) 
is the phase velocity, which is strongly related to f r , 

V = y/I/C (6) 



is the characteristic impedance, L is the inductance per 
unit length, and C is the capacitance per unit length. 
From these transmission-matrix elements, the scattering- 
matrix elements are calculated, and 52 1 is given by 



5 21 = 2/{A + B/Z + CZ + D), 



(7) 



where Zq = 50 £1 is the characteristic impedance of the 
microwave lines connected to the resonator. Unit-length 
properties of CPW are determined when two parameters 
out of v p , Z cpw , L, and C are specified. In the calcu- 
lations for Fig. [4l we employed C = 1.6 x 10~ 10 F/m 
based on the considerations described in the following 
paragraph, and evaluated C c and v p by least-squares fit- 
ting. 

Wer>i£ calculated CPW parameters using conformal 
mapping. Within the theory, C does not depend on t, 
and it is given by 



C= {e r + l)e 2K(k)/K(k'), 



(8) 



where e r is the relative dielectric constant of the sub- 
strate, eo = 8.85 x 10~ 12 F/m is the permittivity of free 
space, K(k) is the complete elliptical integral of the first 
kind, the argument k is given by 



w/(w + 2s) . 



(9) 



and kl = yl — k 2 . For our CPWs, we obtain C = 
1.6 x 10~ 10 F/m when we employ e r = 11.7 for Si (p. 223 
of Ref. d) neglecting the contribution from the Si02 
layer, which is much thinner compared to w, s, or h. Cir- 
cuit simulators [Microwave Office from AWR (#1) and 
AppCAD from Agilent (#2)] also predict similar values 
of C. The simulators calculate CPW parameters from 
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TABLE II: Dependence of coplanar-waveguide parameters on 
the film thickness t. For capacitance C and inductance L 
per unit length, the normalized variations AC(t)/C* and 
AL(t)/L* are listed in percent, where AC(t) = C(t) — C* , 
C ~ 1.6 x 1CT 10 F/m is the value at t = 0.3 fj,m, and the 
definitions of AL(t) and L* ~ 4 x 10~ 7 H/m are similar. The 
predictions by circuit simulators #1 and #2 are compared. 
Regarding L, experimental values for "A" ^Resonators Al- 
A4 and for "B" =Resonators B1-B4 are also given, and they 
are obtained from the values of v p in Table [T] using Eq. ([5]) 
and by neglecting the t dependence of C. 
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the dimensions and the material used for the substrate. 
The predictions by the simulators have t dependence, but 
in the relevant t range, the variations are on the order of 
1% or smaller as summarized in Table [XXI and the values 
of C are between 1.6 x 10" 10 F/m and 1.7 x KT 10 F/m. 
Thus, partly for simplicity, we used C = 1.6 x 10 -10 F/m 
for all of our resonators. 

With C = 1.6 x 1CT 10 F/m, the values of v p in Ta- 
ble [J correspond to Z cpw = 49 — 53 il, which agrees with 
our design of ~ 50 il. We have done the same fitting 
by changing the value of C by ±10% as well in order to 
estimate the uncertainties, which are also listed in Ta- 
ble UJ Within the uncertainties, the values of C c from 
the same coupling-capacitor design agree, and C c ~ 7 fF 
for Resonators A1-A4 with If —78 fj,m and C c ~ 5 fF 
for Resonators B1-B4 with If = 38 /xm. The uncertain- 
ties for v p is much smaller, < 0.1%, and again within the 
uncertainties, the values of v p for the same t agree. 

For the rest of this paper, let us assume that t depen- 
dence of C is negligible. This assumption is consistent 
with the fact that the experimental C c vs. t in Table UJ 
does not show any obvious trend. Moreover, according 
to the circuit simulators in Table UT1 t dependence of C 
is smaller than that of L. Below, we look at L mainly 
instead of v p or other CPW parameters so that we will 
be able to discuss the kinetic inductance. As long as we 
deal with a normalized inductance such as the ratio of 
L(t) to L* = L(0.3 jtxm), what we choose for the value of 
C does not matter very much because v p obtained from 
the fitting was not so sensitive to C. Hence, we ana- 
lyze the quantities obtained with C — 1.6 x 10~ 10 F/m 
only hereafter. In Table UT1 we list the variations of L in 
our two series of resonators as well. For both series, the 
magnitude of the variations are much larger than the pre- 
dictions by circuit simulators. We will discuss this large t 
dependence in terms of kinetic inductance in Sec. IIVI after 
examining the temperature dependence in Sec. IIII BI 
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FIG. 5: (Color online) Amplitude of the transmission coef- 
ficient S21 as a function of frequency for Resonators Al at 
different temperatures. 
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FIG. 6: (Color online) Normalized quality factors, (a) 
Q/Qe and (b) Ql/Qe, as functions of temperature for Res- 
onators A1-A4 (Nb thickness t = 0.05, 0.1, 0.2, and 0.3 (im), 
where Ql, Q e , and Q are loaded, external, and unloaded 
quality factors, respectively, and Q e is assumed to be temper- 
ature independent. The markers are data points, whereas the 
curves are guides to the eyes. 



B. Temperature dependence of S21 

We also measured S21 vs. / at various temperatures 
up to T — 4 — 5 K for Resonators A1-A4. We show 
the results for Resonator Al in Fig. [5] With increas- 
ing temperature, f ri Ql, and the peak height decrease. 
As in Sec. IIII Al let us look at the quality factors first. 
In our resonators, Q e ~ Ql at the base temperatures 
as we pointed out in Sec. IIII Al Thus, when we assume 
that Q e is temperature independent, we can calculate Q 
from measured Ql using Eq. ([T]). We plot QL{T)/Q e 
and Q(T)/Q e vs. T in Fig. [S] for all of the four res- 



5 




12 3 4 
Temperature (K) 



FIG. 7: (Color online) Temperature variations of inductance 
L per unit length for Resonators A1-A4, whose Nb thickness 
is t = 0.05, 0.1, 0.2, and 0.3 /im. The unit of the vertical axis 
is percent. See text for the definition of AL(t,T)/ L* . 



onators. With increasing temperature, Q decreases in 
all resonators. A finite Q^ 1 means that the resonator 
has a finite internal loss, which is consistent with a peak 
height smaller than unity in Fig. [5] The internal loss at 
high temperatures must be due to quasiparticles in the 
superconductor, as discussed in Ref. la. The reduction of 
quality factors becomes larger as the Nb thickness is de- 
creased. At T < 1 K, however, the reduction is negligibly 
small, and thus, in this sense, it should be fine to choose 
any thickness in the range of t = 0.05 — 0.3 /im for the 
study of superconducting qubits that we mentioned in 
Sec. U because qubit operations are almost always done 
at the base temperatures. 

When CPWs are no longer lossless, j3 in Eq. (j4)) has 
to be replaced by (a + 20)1 3- This a characterizes 
the internal loss, and 3 /{2a) is equal to Q (for exam- 
ple, Sec. 7.2 of Ref. 1). From similar calculations to 
those in Sec. IIII A[ we evaluated L at higher tempera- 
tures as well by neglecting the T dependence of C and 
C c . Because we are interested in the temperature vari- 
ation of L, we show AL{t,T)/L* vs. T in Fig. [71 where 
AL{t,T) ee L{t, T) - L(t,T*), T* is the base tempera- 
ture, and L* = L(0.3 (im,T*). The variation becomes 
larger as the Nb thickness is decreased. This trend also 
suggests that we should take into account the kinetic in- 
ductance. 



IV. DISCUSSION 

The film-thickness and temperature dependence that 
we have examined in Sec. IIIII is explained by the model, 



L{t, T) = L g {t) + L k {t, T), 



(10) 



where L g is the usual magnetic inductance per unit 
length determined by the CPW geometry and L k is the 
kinetic inductance of the CPW center conductor per unit 
length. We neglect the contribution of the ground planes 
to L k because the ground planes are much wider than 
the center conductor in our resonators [see Eq. (fTTj) ] . We 
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FIG. 8: (Color online) (a) Temperature dependence of the 
penetration depth A for Resonators A1-A4, whose Nb thick- 
ness is t — 0.05, 0.1, 0.2, and 0.3 (im. The solid and broken 
curves are theoretical predictions expressed by Eqs. (I12|l and 
1(13)1. respectively. As in Ref. 17, [\{t,0) / \{t,T)] 2 is plot- 
ted vs. T/T c {t). (b) Superconducting transition temperature 
T c {t) used in (a), and \{t, 0). The curves are from Ref. Il2l. 
that is, not fitted to our experimental data. Both in (a) and 
(b), \{T*) ~ A(0) is assumed, where T* is the base tempera- 
ture. 



also assume that L g depends on t only, whereas L k does 
on both t and T. This type of model has been employed 
in earlier work o 9 i 13 ' 14 as well. The T dependence of L k 
arises from the fact that L k is determined not only by 
the geometry but also by the penetration depth A, which 
varies with T. Meservey and Tedrow 8 calculated L k of a 
superconducting strip, and when the strip has a rectan- 
gular cross section like our CPWs, L k is written as 



^r{\/w) ln{4w/t) 



sinh(i/A) 
cosh(t/A) - 1 



(11) 



where fiQ — Air x 10~ 7 H/m is the permeability of free 
space. The relationship between Lk and A is expressed 
in a much simpler form in the thick- and thin-film limits; 
L k cx A for t ^> A, and L k cx A 2 for t -C A. When we 
assume Eqs. (ITU)) and (fTTj) . we obtain A(t, T) numerically, 
once L g {t) is given. Below, we discuss \{t,T) in our Nb 
films in order to confirm that the model represented by 
Eq. (|10|) is indeed appropriate. 

In Fig.HHa), we plot [A(t, T*)/X{t, T)] 2 vs. T/T c {t) for 
Resonators A1-A4, where T c {t) is the superconducting 
transition temperature, which is assumed to be also t 
dependent in this paper. We have found that with a 
reasonable set of parameters, L g {t) and T c {t), the exper- 
imental data for all resonators are described by a single 
curve. This kind of scaling is expected theoretically in 
the limits of £o/Al S> 1 and £oMl <C 1, where £o is 
the coherence length and Ai is the London penetration 
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TABLE III: Inductance per unit length at the base temper- 
atures in Resonators A1-A4. t is the thickness of Nb film; 
AL g (t) — L g (t) — Lg, where L g is the usual magnetic in- 
ductance per unit length determined by the CPW geometry, 
and L* g = L 9 (0.3 (im) = 3.75 x 10~ 7 H/m; Lk is the kinetic 
inductance per unit length, and L = L g + L^. 
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depth Although £, /X L ~ 1 in Nb (p. 353 of Ref. QJ) 
at temperatures well below T c , it would be still reason- 
able to expect a scaling in our Nb resonators because 
at a given normalized temperature T/T c (t), the relevant 
quantities should be on the same order of magnitude in 
all resonators, and thus, two parameters, A(t, T*) and 
T c (t), are probably enough for characterizing X(t,T) of 
our resonators. The values of L g (t) and T c (t) employed 
in Fig. [UJa) are summarized in Table IIIII and Fig. [8jb), 
respectively. The relative change of L g (t) in Table IIIII 
is similar to the predictions by circuit simulators in Ta- 
ble HH which do not take into account the kinetic induc- 
tance. The magnitude of L g (t) is also reasonable because 
y/L g (t)/C ~ 49 fl for all thickness. In Table [ml we 
also list the ratio of kinetic inductance Lk to the total 
inductance L. With decreasing thickness, Lk/L indeed 
increases rapidly. In Fig. IDJb), T c (t) and X(t, T*) are 
plotted together with the theoretical curves in Figs. 1 
and 6 of Ref. [ID, where Gubin et alw^ determined some 
parameters of the curves by fitting to their experimental 
data. The values of T c (t) are reasonable, and X(t,T*) is 
on the right order of magnitude. 

The solid curve in Fig. [Sf a) is the theoretical T depen- 
dence based on the two-fluid approximation,— 

[A(0)/A(T)] 2 = 1 - (T/T c )\ (12) 

This theoretical curve reproduces the experimental data 
at T/T c < 0.4, when we assume that X(t, T*) — A(i, 0) in 
Resonators A1-A4. At T/T c > 0.4, on the other hand, 
the experimental data deviate from Eq. (fT2")l , but accord- 
ing to Ref. Qj], the expression for A vs. T depends on the 



ratio of ^o/Al, and thus, Eq. (|12p cannot be expected to 
apply to all materials equally well. Indeed, although the 
temperature dependence of Eq. lfl~2"j) has been observed 
in the classic pure superconductors?^ such as Al with 
Co/A/, ^> 1 at temperatures well below T c , it does not 
seem to be the case in the high-T c materials, whose typ- 
ical £o/Al is in the opposite limit ^ £o/Al <§; 1, and for 
example, Rauch et al. 9 employed for a high-T c material 
YBa2Cu307_j;, an empirical expression of 

[A(0)/A(T)] 2 = 1 - 0.1(T/T C ) - 0.9(T/T C ) 2 , (13) 

which is the broken curve in Fig.[5Ja), instead. Because 
Co/Al ~ 1 in Nb even at T/T c <C 1, and because the 
experimental data at T/T c > 0.4 are between Eqs. (fT2)) 
and (fTB")) . we believe that the deviation from Eq. fT2")) at 
T/T c > 0.4 is reasonable. 

From the discussion in this section, we conclude that 
the model represented by Eq. (|10p explains the film- 
thickness and temperature dependence of our resonators. 



V. CONCLUSION 

We investigated two series of Nb A/2 CPW resonators 
with resonant frequencies in the range of 10 — 11 GHz 
and with different Nb-film thicknesses, 0.05 — 0.3 /xm. 
We measured the transmission coefficient S21 as a func- 
tion of frequency at low temperatures, T — 0.02 — 5 K. 
For each film thickness, we determined the phase veloc- 
ity in the CPW with an accuracy better than 0.1% by 
least-squares fitting of a theoretical £21 curve based on 
the transmission matrix to the experimental data at the 
base temperatures. Not only the film-thickness depen- 
dence but also the temperature dependence of the res- 
onators are explained by taking into account the kinetic 
inductance of the CPW center conductor. 
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